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The techniques used previously to calculate the three-fold thermal conductivity reduction due to 
phonon dispersion in GaAs/AlAs superlattices (SLs) are applied to HgTe/CdTe SLs. The reduction 
factor is approximately the same, indicating that this SL may be applicable both as a photodetector 
and a thermoelectric cooler. 
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I. INTRODUCTION 

Hg^Cdi-^Te alloys are among the key materials used 
for cooled infra-red detectors. Materials structures in- 
volving these constituents in both alloy and supcrlat- 
tice (SL) forms that would be useful for both opti- 
cal detection and thermoelectric cooling to a reason- 
ably low temperature are an attractive possibility. Some 
semiconducting superlattices (SLs) are known to exhibit 
an order-of-m^i^u£^f«d-action in the lattice thermal 
conductivitie^i'BMmGl'BBM. For example, Capinksi et 
al have experimentally shown GaAs/AlAs SLs to have 
a ten-fold smaller lattice thermal conductivity than the 
bulk alloys having the same compositiorOllj. This result 
implies that the thermoelectric figure of merit ZT would 
be increased by an appreciable factor. We show theoret- 
ically that the same appears to be true for HgTe/CdTe 
SLs. Cooling of integrated room temperature devices to 
200-250 K thus becomes a possibility for applications re- 
quiring no external coolant provided these temperatures 
are adequate. Since the devices would be fabricated 
of the same material constituents, manufacturing costs 
might be substantially reduced. 

We have-, previously investigated GaAs/AlAs 
theoreticalliJip using Kunc's model for the lattice 
spectrumEJllj generalized to SLsEj. Details are given 
in Ref. 13. This rigid-ion model considers next-nearest 
neighbor short-range forces together with the long-range 
Coulomb force. We found a factor-of-three reduction 
in the lattice thermal conductivity along the growth 
direction due to SL induced modifications of the phonon 
spectrum. The remainder of the overall ten-fold reduc- 
tion is presumably associated with phonon scattering 
at the interfacesEJ'OEj which is not considered in the 
present work. The thermal conductivity reduction 
in question here can be interpreted as resulting from 
zone folding in the SL relative to bulk. The new gaps 
introduced at the zone face result in band flattening 
and thus lower phonon velocities, which in turn reduce 
the thermal conduction. These effects are illustrated in 
Fig. 1 and Fig. 4 of Ref. 13. The corresponding figures 
in this paper will be seen to exhibit very similar effects 



in HgTe/CdTe SLs. For reasons to be discussed, these 
results lead us to expect similar behavior in any SL 
having the same structure and similar mass and force 
constant differences between layers. Due to residual 
"impedance mismatches" at the interfaces, the Kunc 
model leads to decaying TA and LA modes along some 
intervals and directions in q-space. (This is not the 
case for GaAs/AlAs SLs.) We resolve this problem by 
use of a variety of physically motivated modifications, 
which are critically examined in the next section. The 
same three-fold reduction previously calculated for 
GaAs/AlAs is obtained. Moreover the result is shown 
to be robust with respect to the types of physically 
reasonable approximations made. Our discussion em- 
phasizes that a reasonable phonon spectrum is obtained 
if one of the force constants, obtained from neutron 
diffraction experiments for HgTe and CdTe respectively, 
is modified by a small amount. Unfortunately there are 
no corresponding neutron data for the SL that might 
lead to better experimental input for determining the 
Kunc force constants. 



II. THE KUNC MODEL FOR HgTe/CdTe SLs 



The Kunc rigid-ion modeli3't3 has been successfully 
applied to compounds with zinc-blende structure. In 
each bulk compound, ten parameters describe the force 
constants between nearest and next-nearest neighbors 
and one parameter the effective charge. In the absence 
of experimental data for HgTe/CdTe SLs there are no 
directly determinable parameters with which to fit the 
phonon spectrum. Note that even for bulk HgTe the 
fit to the neutron-diffraction data along the [111] direa- 
tion does not reproduce the optical branches very wellta. 
Thus, we use the parameter fits to bulkpHgTe and CdTe 
determined by Talwar and VandevyveitJ and construct 
the SL dynamical matrix from these. There are thus 22 
parameters needed for the SL, in addition to the atomic 
masses and lattice spacing. The bulk parameters are used 
within each layer. At the interface we average the force 
constants from the two sides. The Madelung sum is cal- 
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culated using the Ewald transformation adapted to the 
SL. We refer to Ref. [l3] for the details. The output of 

the Kunc model is the phonon dispersion o; q Q ^ for the 
a-th branch as a function of wavevector q in the SL Bril- 
louin zone. The lattice thermal conductivity may then 
be computed from the phonon Boltzmann equation in 
the relaxation-time approximation: 
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where n(oj q Q ^) is the distribution function of phonons and 

Tph(w q a ' ) , T) is the lifetime, taken to be a constant r. 

The dynamical matrix for the HgTe/CdTe SL so con- 
structed, though Hermitian, has negative eigenvalues, 
i.e., yields imaginary phonon frequencies, for the three 
acoustic branches in a cylindrical region of the Brillouin 
zone around the [001] axis that occupies about 6% of the 
zone by volume. The problem of imaginary frequencies in 
models of HgTe/CdTe SLs has been investigated by Ra- 
jput and BrowneEj They attribute the problem to an in- 
adequate modelling of the interface, in their simple force 
model based on the Keating potential. They find that 
in their model the imaginary frequencies can be removed 
by treating the two materials as continuum media with 
differing dielectric constants. However, this prescription 
fails to remove the imaginary frequencies using the Kunc 
model. As we shall see below, the imaginary frequencies 
in the Kunc model are more likely to arise from the force 
constants than from the Coulomb interaction. 

In order to gain an idea of what is causing the imag- 
inary frequencies, we looked at the dependence of the 
gap at the zone edge, Aw, on the Kunc parameters. We 
found the strongest dependence for the diagonal terms in 
the force constant matrices; the effect of the off-diagonal 
terms was weaker by a factor of ten or more. This finding 
motivated us to consider the variation of the parameters 
necessary to remove imaginary frequencies. For example, 
an interpolation of all parameters (20 force constants, 2 
effective charges, 3 atomic masses and the lattice spacing) 
linearly between between their values for a HgTe/CdTe 
SL (x — 0) and those for a GaAs/AlAs (x — 1), which 
has real frequencies only, shows that all of the imag- 
inary frequencies disappear when x approaches 0.035. 
This suggests that we associate the imaginary frequen- 
cies with an impedance mismatch between the HgTe and 
CdTe layers which can be removed by variation of appro- 
priate force constants or the introduction of additional 
ones provided by a more sophisticated model. Thus, if 
we hold the parameters fixed at their HgTe/CdTe values 
and vary the parameters individually, we find that small 
(< 10%) changes in An g Te and AcdTe alone are suffi- 
cient to yield real frequencies throughout the SL Brillouin 
zone. Here, A is the nearest-neighbor coupling between 
Hg and Te respectively Cd and Te. Apparently the prob- 
lem does not exist in GaAs/AlAs SLs because the A pa- 



rameters are so closely matched: AcaAs = —40.77 N/m 
and Aaias = -40.83 N/m. 

We now describe our approaches for removing the 
imaginary frequencies. The only correct way to solve the 
problem is to generalize the already complicated Kunc 
model. However, the Kunc-model description of the bulk 
phonon dispersion is satisfactory. Changing the bulk 
model does not necessarily imply that the SL results will 
be improved because the interfaces are neglected. In- 
stead, we construct physical models, for example, deter- 
mining values of A that eliminate the imaginary frequen- 
cies. The reduction factor is always about 3 for the more 
realistic models. The robustness of the result gives us 
more confidence in its physical relevance than would de- 
riving it from a generalized Kunc model with many new 
parameters. 

The best model utilizes scaling, which maintains the 
acoustical and optical mode frequencies observed in bulk. 
The effect of changing A, aside from removing the imag- 
inary frequencies, is to rescale all the phonon frequencies 
across the Brillouin zone. The scaling model is therefore 
constructed as follows: (1) A is replaced by A'; in prac- 
tice A 1 is selected so that the ratio of LA to TA sound 
velocities in the [001] direction will be that of bulk; this 
means x — 0.2 or AngTe is changed from —18.127 N/m to 
-22.7 N/m and A C dT C from -19.97 N/m to -24.1 N/m. 
This gives sound velocities about 10% too high and in- 
creases the frequencies of the optical branches by 10%. 
(2) In order to compensate for this increase, all phonon 
frequencies are rescaled by a common factor, lu — ► Xlu, 
with A = 0.927, assuming A to be constant. The scal- 
ing model has the advantage of staying within the Kunc- 
model framework, while removing the imaginary frequen- 
cies in a consistent way. 

Another model is based on averaging the force con- 
stants (both nearest neighbor and next-nearest neighbor) 
at the interface. This approach serves to equalize the 
force experienced in each direction by the interfacial Te 
ion. This approach fails, however, to remove the imagi- 
nary frequencies. The fraction of the Brillouin zone cov- 
ered by imaginary frequencies changes by less than 1%. 

The problem can be addressed most simply by setting 
the frequencies of the imaginary branches equal to zero. 
We find a thermal conductivity reduction by a factor of 
3.0, very close to that found for GaAs/AlAs. Since only 
three of the eighteen folded-back acoustic branches are 
set to zero over a small part of the zone, the end re- 
sult should be close to the true value, which represents 
a lower bound to the reduction factor. For comparison, 
we find that setting the lowest three branches to zero in 
GaAs/AlAs gives a reduction by a factor of 2.9, versus 
2.8 in the correct calculation. Thus we can expect that 
the actual value in HgTe/CdTe will be close to 3. 

Another approach involves a linear interpolation from 
the origin to the point of tangency with the acoustic 
branches where they are real. This gives a reduction 
factor of 2.9. The treatment of the acoustic branches 
however is unrealistic because the interpolated sound ve- 
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locities are too low, only about 50% of the bulk value 
(average of HgTe and CdTe) in the [001] direction. 

The simplest model is based on the virtual crystal ap- 
proximation, in which the force constants are replaced 
with their averages between the two layers, while the 
masses and effective charges are kept at their respective 
values for the two layers. This ensures that the frequen- 
cies are real, because of the absence of interfaces. How- 
ever, it does not properly account for the thermal conduc- 
tivity reduction (giving only a factor of 1.3), because the 
band gaps at the zone edge, which reduce the phonon 
velocity, are considerably smaller than they should be. 
These results indicate that it is the force constants and 
not the Coulomb interaction that are responsible for the 
imaginary frequencies in the SL Kunc model. 



III. RESULTS 

We consider a (HgTe)3/(CdTe)3 SL. The dispersion re- 
lation along the FX and TZ directions was generated nu- 
merically using the scaling model as described in Sec. II 
and is shown in Fig. 1. The SL unit cell contains three 
unit cells each of HgTe and CdTe, arranged along the 
growth axis. Thus, the edge of the SL Brillouin zone 
is one-sixth as far from the center along the TZ axis 
as it is in the in-plane TX and TY directions. Each of 
the six branches in bulk is folded back six times in the 
TZ direction, yielding 36 branches in the SL. This ef- 
fect is most easily seen for the acoustic modes along TZ. 
The results are qualitatively similar to those obtained 
for (GaAs) 3 /(AlAs) 3 in Ref. [l| The HgTe and CdTe 
optical modes are not as well separated as the GaAs and 
AlAs optical modes are in GaAs/ Al As, and we see that 
the lowest HgTe optical modes overlap with the high- 
est acoustic modes along FX. Nevertheless, the optical 
modes are nearly flat in the TZ direction. 

As can be seen in Fig. 1, the physical effects of band 
flattening, relative to bulk, are present. We therefore 
expect, by Eq. ([l]), a reduction in the thermal conduc- 
tivity along the growth axis. In fact, with our data, 
Eq. (|l|) yields a reduction factor in k/t of 3.0 (in the 
scaling model). This factor reflects phonon-dispersion 
effects only. As discussed in Ref. [l3|, one may expect a 
further reduction in thermal conductivity if the lifetime 
r itself is lowered in the SL compared to bulk due pri- 
marily to interface scattering. But there is not sufficient 
data on HgTe/CdTe SLs for us to discuss the lifetime 
quantitatively. 

A simple physical picture of the thermal con- 
ductivity reduction due to phonon dispersion effects 
is obtained by looking at the transport quantity 

<?z EJ*(4°')/ (ir ) u q a) (''v) 2 in q space. This is the 
quantity in square brackets in Eq. ([!]). Since the dis- 
persion relation is rotationally symmetric in the (q x ,q y ) 
plane to good approximation we perform an annular inte- 
gration, yielding the factor q x equal to the annular radius 
and leaving an integral in Eq. ([!]) depending only on q x 



and q z . The resulting quantity is plotted in Figs. 2(a), 
2(b) and 2(c) for the SL and for bulk HgTe as a function 
of q z for three values of q x . The SL transport quantity is 
seen to be localized, i.e., it tends to zero at the band edges 
q z = and q z = it / d. The reduction in the SL transport 
quantity near the zone edges is related to miniband for- 
mation (the flattening of a; versus q z seen in Fig. 1). To 
gain physical insight, we replace the SL transport quan- 
tity in Figs. 2(a), (b) and (c) by rectangles of equal area 
and plot these in Fig. 2(d), which summarizes the depen- 
dence of the transport quantity on q x . The shading indi- 
cates the weight of each equivalent rectangle within the 
range {q Xl q x + Ac^). In bulk the equivalent rectangles 
extend all the way to the zone edges because there is no 
localization due to band flattening. In Fig. 2(d) we see 
that the densities tend to zero as q x goes to zero and 
that the the density in bulk increases for large q x , again 
due to the q x prefactor. Surprisingly, the density in the 
SL falls off fast enough as q x increases to overcome the 
q x factor. This should be compared with the results for 
GaAs/ AlAs SLs of Ref. The qualitative difference be- 
tween HgTe/CdTe and GaAs/ AlAs SLs is not apparent 
from their dispersion relations. Integrating the density in 
Fig. 2 gives an estimated reduction factor of 2.4, which 
is to be compared with 3.0 in the exact calculation. The 
difference is due to the fact that the estimate does not 
properly treat the region of the Brillioun zone with ra- 
dius greater than \Z2ir/ao and therefore underestimates 
the reduction factor. 

Our results predict a reduction factor in ther- 
mal conductivity for HgTe/CdTe SLs comparable in 
size to those already observed for GaAs/ AlAs SLs. 
The minimum cold temperature achievable for a 
30A Hgo.75Cdo.25Te/30A Hgo.7Cdo.3Te SL is about 180 
K, compared to 240 K for conventional Bi2Te3. This 
result suggests that it is reasonable to use HgCdTe for 
thermoelectric cooling applications in infra-red detectors. 
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FIG. 1: (HgTe)a/(CdTe)3 superlattice dispersion relation 
along the rX=(27r/ao, 0, 0) and TZ=(0, 0, 7r/3ao) directions, 
where do = 6.48 A is the conventional unit cell size of bulk 
HgTe. 
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FIG. 2: (a),(b),(c) The transport quantity 

{dn/dT)uq"\vqfl) 2 for < q z < n/d at fixed q x 
for (a) Qa; = 2n/d, (b) = 4n/d and (c) ^ = 6ir/d. The 
bulk HgTe and the (HgTe) 3 /(CdTe) 3 SL curves are labelled 
on the right of figures (a), (b), and (c). Note the smallness of 
the SL curves in (b) and (c). (d) Density plot in the {q x ,q z ) 
plane whose shading indicates the weight of each increment 
in Aq x along q x to the value of the transport quantity for the 
(HgTe) 3 /(CdTe) 3 SL and bulk HgTe respectively. 



